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STUDY OP A PRECISE POSITIONING SYSTEM 
OP THE COMSS ALTIMETRIC SATELLITE/RADAR 
(HASP) 

(High Accuracy System Positioning) 


P. Nouel? M. Berge** 


1. OBJECTIVES OP THE STUDY /6«** 

If a satellite is equipped with an altimeter whose measure- 
ments have an accuracy on the order of between E to 10 centimeters 
(a noise), it is apparent that the satellite orbit must be 
calculated with an Identical accuracy level or better, and this 
is practically indispensable in the case of the SEAST. In order 
to solve calibration problems of the altimeter, the exactness 
with which the trajectory is reconstructed also is Involved. 

Therefore the proposed system from this study will respond 
to the first objective: reconstruction of the satellite tra- 
jectory at an accuracy level of E-lO centimeters. 

For this purpose, two other objectives have to be taken 
into account: possibility of reconstruction of the orbit at a 
level of 1 to ? meters, but during a shorter time Interval, 
l.e., between 1 to ? days; comparison of the entire results for 
two different altitudes, 6E0 and SEO km. 

f. HOW CAN WK PROCEED WITH THE STUDY? 

The first difficulty which one encounters is to know how 
such a study should be carried out. In effect, how can one 
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estimate 5 centimeters (Figure 1) In the position of a satellite 
which is at an altitude of 600 km or more, and which traverses 
a distance of more than 40,000 km during one orbit around the 
Earth? 

On the other hand, the HASP system does not exist, and 
thus no constraint has to be Imposed on the methods considered. 
The study will be based on tracking techniques which now exist, 
or which probably will exist very shortly, and this will be 
associated with processing techniques. This association has 
rarely been taken Into account In the past. 

Considering the experience acquired and the enormous effort 
expended for calculating the SEASAT trajectory, the directions 
of research are relatively well defined. 

It is within this context that the study was carried out. 
The uniqueness of the system will depend only on the initial 
choice of the research method. 

3. GENERAL REMARKS ABOUT THE TRAJECTORY CALCULATION 
3.1. Vector representation 

No matter what type of measurements are considered, the 
satellite station vector ( q ) is Involved (the two Instruments 
at each extreme point are assumed to be points) (Figure 2): 

. either in terms of direction (photograph agalns^ a star 
background, interferometry), 

. or by its modulus (distance measurement) 

. or variations in its modulus (Doppler measurements) 


This vector will be calculated in an inertial reference /8 


T^lgur« 1. 5 centimeters! 



(measured ) 
ground station 


mass center of 
the earth 


Figure 2 


i 


3 


i 


I 


system in order to write the equations of dynamics. 

Therefore, at the station extremity it is necessary to know 
the position of the station (15 in a coordinate system, and 
this system has to be referred to an Inertial system, 
obvious that the accuracy with which one can transfer 
system to another will evolve from this study. 

As far as the satellite extremity is concerned 
motion is due to a certain number of forces for which 
exist, but which are more or less representative. As 
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example, we have the following: 


- from the attraction due to the Earth’s potential. The 
Earth's potential Is usually extended In a series of 
spherical harmonics, 

- gravitational attraction forces from the Moon and the 
Sun which are sufficiently well known for satellites 
close to the Earth, 

- attraction forces caused by the potential of terrestrial, 
oceanic, and atmospheric tides, which Introduces long 
term perturbations, 

- the force of atmospheric drag which obviously Is limited 
to a certain altitude range, 

- forces caused by pressure and direct radiation of the Sun 
and the pressure re-emltted by the Earth, In particular. 

In the case of a hello-synchronous orbit, the direct 
pressure will have a relatively constant and continuous 
effect. The re-emltted pressure becomes variable for 
shadow-Sun passages, and for cloud covers. This is 
difficult to model. 

3.2. Calculation techniques 

At the present time, calculation techniques for following a 
satellite consist of calculating an orbit which represents the 
measurements the best. This orbit Is then used to interpolate 
or predict the satellite positions. Considering the perturbations 
and the required accuracy, this evaluation may not be done 
to a great extent: this is because the trajectory arc lasts for 
several days. 

In the following study, we shall use the notation below: 

short arc - 

Thls Is one part of the satellite orbit which can extend 
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up to one complete revolution around the Earth. 


Long arc ; 

If we consider several successive revolutions and even 
orbits which extend one or several days, we will be dealing with 
a long arc. Considering the desired accuracy for space 
geodesy for a close satellite. In general the duration will not 
exceed 5 days. 

EVALUATION OP THE FORCES ACTINO ON A SATELLITE 
FOR A SHORT ARC 

Approximately, a close satellite of the Earth traverses 

7.000 km In ?0 minutes (about 1,000 seconds). Thus an ac- 

-7 2 

celeratlon of 10 m/s will result In a displacement on the 

order of 5 centimeters In 20 minutes, which is within the 

range of the present study. Consequently, the accelerations 

from which the satellite motion will be calculated must be known 
-8 2 

to 10 m/s^*, so that the error remains smaller than 1 centimeter 

over a period of 20 minutes. /lO 

Non-gravl tat lo nal forces 

In order to order the orders of magnitudes as a function of 
altitude for atmospheric friction, for the same satellite 
(example taken from D5B), we will have an acceleration (in 
m/s^") of 10”^ at 300 km, lO”"^ at 600 km, and 5*10“^^ at 

1.000 km. 


The following are accelerations due to pressure forces: 

- 3'10~ for direct solar pressure 

-9 

- 3’ 10 for reflected solar pressure 

- ^•^0”'^ for the Inrared. 


5 


The problem arises because one encounters variations by a 
factor of 100 of these accelerations during an orbit (passage 
into the shadow), and even more so over a day. 

Forces caused by terrestrial potential 

When a new geodesic satellite Is launched, the first 
factor encountered when one determines the orbit Is that the 
potential use Is Inadequate. 

For example, we can mention GEOS III which led to determina- 
tion of the models GEM 9 and 10. When SEASAT was launched, it 
was necessary to calculate the models PGS-Sl and PGS-S2, and 
we not only did not have an improvement in the arcs from one 
to three days, but also for arcs lasting several hours. 
Nevertheless, the geographic distribution of observation stations 
is an important factor in determining most of the coefficients 
of the potential model. For various reasons (isolation, political) 
certain regions of the globe will remain inaccessible. /II 

5. PRESENT ACCURACY OP ORBIT CALCULATIONS: SEASAT EXAMPLE 

Prom known results for several geodesic satellites, l.e., 
where it was desired to obtain the best accuracy in the calcula- 
tions, we can make the following estimates: 

The radial component Is calculated with an error between 
1 to 2 meters In general. If particular efforts are expended, 
we can reach an accuracy of 50 centimeters. 

This accuracy is valid for arcs varying between 1 to 2 
revolutions up to 1 to 2 days In the more precise case. 

The SEASAT results In addition show the following predominant 
error sources due to our lack of knowledge: 
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- terrestrial 

potential 

30?S 

- atmospheric 

friction 

20 % 

- solar pressure 

20 % 

- propagation 

phenome.ia 

10 % 

- station positions 

10? 

- calculation 

procedures 

10? 


6. METHOD OP "SHORT ARCS" 

In accordance with a goal of space geodesy, new methods of 
calculation have been proposed In order to calculate the 
positions relative to the stations on the ground. 

The general philosophy of these methods consists of 
observing the same satellite trajectory arc from several stations 
so that the errors Introduced by the poor knowledge of the 
satellite position are common to all of the stations, and thus 
their effects are re-^'iced. /\ 

For example, with the satellites TRANSIT, whose accuracy 
of the trajectory Is on the order of 25 meters, If one uses 
ephermerldes which are called operational and which are 
transmitted by the satellite, we can calculate the station 
positions with an accuracy on the order of 50 centimeters to 
1 meter. 


Another very characteristic example is the estimation 
of the quality of measurement apparatus: the apparatus noise 
can be estimated by using measurements collected over a 
pass above the observation station (case of laser with an 
accuracy between 2 and 5 centimeters) with data processing 
whose satellite motion model is simple. 

This method of short arcs has one drawback: this increases 
the number of measurement apparatus, but the measurement 
techniques have Improved accuracy, and as a result long wave 
radio waves have to be used. 


7. EXISTENCE OF TRACKING MEANS 


At the present time, laser techniques give the required 
accuracy (2 cm) in distance, but in spite of the efforts to 
reduce their size (mobile stations and very mobile stations), 
the construction, installation, and operation have Increased 
payoffs . 

On the other hand, we have seen a very Important improve- 
ment in radio techniques, for example: 

- the Transit system: Doppler receivers now have a small 
size 

- automation of measurements: Eole or Argos. /13 

- appearance of distance measurements with an accuracy and 
exactness of several centimeters (Dialogue model) or 
systems with real time exploitation (GPS system). 

- VLBI stations which are mobile and are still in the 
experimental stage and are still quite large 

The operational problems and cost problems have partially 
been solved, and at the present we can look upon tracking equip- 
ment in a different light. 

A network which is more dense could be considered with quasi- 
automatic collection of data and centralization of calculations. 


1. SELECTION OF MEASUREMENT INSTRUMENT 
Present system s 

Laser 


/15 


The measurement accuracy is on the order of 10 to 20 cm, 
and for certain stations it is only a few centimeters, for any 
measured distance. 

The calibration problems, i.e., the measurement of the sys- 
tematic effects which affect the measurements, are relatively 
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easier to resolve than for the radio- location s. Neverthe- 

less, on a centimeter level the laser method encounters the 
same problem. There Is the delay due to the electronics, the 
delay which varies with temperature, the sensitivity of the 
photomultiplier, etc... 

Prom the satellite point of view, the laser requires several 
on-board reflectors, which have minimum cost. On the other hand, 
the ground stations are substantial, even for the new generation 
of mobile stations. 

Finally, their operating cost Is very high. Thus they have 
been considered In numerous programs. 

NASA (Including SAO), Interkosmos, and Europe are almost the 
only organizations which have them. Only NASA and Interkosmos 
have deployed a global network. This Implies international co- 
operation which is difficult to Implement. 

There is still the problem of operations during cloud cover 
which over the year amounts to almost 30% success rate for the 
number of programmed passes: this number varies from one 
geographic location to another. /16 

RADIO LOCATION 

Transit System 

This is based on measuring the Doppler effect alone, and 
is used everywhere throughout the world. It was developed in 
i960, and therefore a long time ago from a technological point 
of view, in spite of important Improvements, especially concerning 
the oscillator stability. It has been successful especially 
in space geodesy, and has been reliable. The costs of operation 
are relatively low, and operating it is relatively simple. 








Obviously, the satellite segment is more complex than in the 
case of the laser. 

This system operates practically for all tines and the 
success rate is between 70 to 80!J when one uses the automatic 
mode, l.e., when all passes are observed, even those having a 
maximum elevation less than 10°, which Is included in 20 to 
30? of rejects. 

OPS System (Global Positioning System) 

It is being installed at the present time with 6 satellites. 

The GPS will give access to distance and Doppler measurements. 

For reliability reasons, the USAF has turned its attention 
towards radio-locating techniques. For the Transit system, the 
most important ^eaturer which were adapted for defining the 
GPS can be summarized as follows: 

- use of a dlstance-Doppler pair which gives a better 
geometric definition of the measured point 

- visibility of several satellites from one point on the /17 
Earth so as to incrase the geometric factor in 

localizing the point \ 

I 

- utilization of higher frequencies in order to reduce the | 

influence of the ionosphere 

- stability of oscillators in order to increase the measure- 
ment accuracy | 

EOLE, ARGOS Systems I 

I 

I 

Here, the desired localization accuracy is not comparable 
with that of the laser or the GPS. This was not the principal 
objective. In the two cases, it was desired to automate the 
localization of s large number of objects (balloons, buoys, etc.) 
which are distributed over the globe. It was desired to 
rapidly give the users the results of the position determinations 




by using the satellite Itself as a data transmission vehicle. 

Even though the concept Is relatively old, it we not widely 
accepted by users which wanted to be autonomous, often a 
fictitious factor. Information processing and telemetry has 
ohangedthls point of view more and more. Centralized data 
bases have appeared, which are acceptable to many users. 

Interferometer systems 

These are based on the phase comparison of signals coming 
from the same source Lud received at two different points. This 
technique allows one to determine the orientation of the satellite 
(DIANE system) . The distance between the stations on the ground 
has been determined better, because of more precise oscillators 
and a better synchronization of the local clocks. Phase /l8 

comparison is carried out at different times on computers. 

Let us recall that, by increasing the distance of the base, we 
can obtain a better angular resolution. Radio astronomers have 
used this method. 

Such a system Is now being developed, and tested, but the 
station Is still costly and substantial, at least in its initial 
design . 

Radio Interferometry as a means of loc/.lizing close satellites 
will develop In the near future, but within a different context 
from the present time. Comparing the phase of a signal to 
anothei phase, of course, leads to a measurement of distance, 
Doppler measurements, and Interferometry. 

CONCEPT OP THE HASP MEASUREMENT INSTRUMENT 

Definition 


For reliability reasons, accessibility, cost, we will not 
discuss lasers. We wish to develop a new instrument without 
being compromised by previous commitments (such as NASA with 
the laser) . We will have to draw on previous experience of radio 


tracking, and use the lastest techniques available In electronics 
and calculation methods. 

It Is always possible to place several laser reflectors on 
the satellite, In order to enlarge International participation 
and In order to verify the accuracy of the proposed sytem using 
another technique. /19 

Orbit coverage 

It Is Imperative to have the best possible orbit coverage: 
example, use of "S-band” measurements for calculating the 
SEASAT trajectory. 

This leads to a radio tracking system. 

Dlstance-Doppler pair 

The measurements of distance and Doppler measurements com- 
plement each other perfectly; the position of the satellite Is 
defined by six orbital parameters, or in an equivalent manner 
by its position and velocity. 

The measurement of distance is performed along the statlon- 
salelllte direction, and the Doppler measurement (or radial 
velocity) Is done In the direction along the line perpendicular 
to the distance. 

Example: Identlfical selection In GPS 

Measurement of double trajectory distance 

In GPS, distance Is defined along a simple trajectory: this 
Is only for military reasons, because the ground station has to 
remain quiet. Consequently, the satellite and the station have 
to both have frequency references and very stable time bases, 
which are set Into phase before they are used. 


12 


In a simple trajectory D (measured) ■ d (geometric) 

♦ c At where c is the speed of light and At is the 
shift between the two clocks. This correction term disappears 
.or a double trajectory. 

Double trajectory Doppler measurement /20 

This problem is simV' to the previous one. TRANSIT or 

GPS carry out simple traj ;ry Doppler measurements. 

One could believe that, in contrast to the distance which 
gives a measurement scale related to frequency, the effect of 
a Doppler frequency deviation does not have any effect. This 
is true if one considers a single measurement, but as soon as 
one wants to mix measurements from several receivers it is 
necessary to know that the frequency deviations are known. 

Simple Doppler 

Double Doppler 

where f^^ is the measured Doppler frequency 

f^ is the receiver frequency 

f is the satellite frequency 
c 

v is the speed of the satellite relative to the station 
c is the speed of light 

Shipping of the measurement apparatus 

The double trajectory principle in distance and Doppler 
can be conceived along the station-satellite-station line: 
the satellite segment will be simplified. 

In this ruse on the ground we multiply the number of 
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reference oscillators, and It is necessary to provide synchro- 
nization (on the order of one microsecond) between all clocks. 
This is the last problem, the problem of transmitting measure- 
ments to the calculation center. This problem is costly, and 
will result in substantial delays in the processing center. /21 

Shipping of the measurement apparatus implies miniaturiza- 
tion in terms of weight and volume, and reduced power consump- 
tion. On the other hai.d, for all of the distance and Doppler 
measurements we have a single frequency reference, and a single 
time reference. This eliminates an Important source of errors 
which are present in present systems. (Even for lasers, where 
Incorrect dating leads to rejection of measurements which 
often have an accuracy of several centimeters!, and this 
comes about because of poor synchronization). 

The last, and not least Important, advantage Is the rapid 
centralization of measurements. 

The groud equipment is simplified, and the number of 
stations can be Increased. 

Summarizing: the following diagram shows the system 
which could be used in HASP: 

. Doppler and distance double trajectory measurements 

. shipping of measurement apparatus 


1 ^< 





1 single clock 
1 oscillator 
Back-and -forth 
distance and 
Doppler 


Measurement transmission 


Calculation center 


5^ ; ^ * Reometrlc satellite positions 

6 A A : beacons 

a, ®a. / 


COFiRECTIONS TO RE MADE TO MEASUREMENTS 

In effect, the station-satellite radius vector 
goes from a station point which is geometrically well defined 
to the satellite center of mass, to which the dynamic equations 
will be applied. 

For each measurement it is necessary to correct either 
geometric effects, or propagation effects wi';h an accuracy 
greater than the !:;-10 centimeters which one desires to re- 
construct . 






SATELLITE 



Ihe list of corrections which we mentioned is not exhaustive 
and we cannot give a solution to all of these problems here. /23 
Each one has to be dealt with by a certain specialist. This 
list Is only Intended to provide the corrections which can be 
predicted and which have to be studied. Certain responses have 
already been made during development. 

Geometric corrections 


- Studies on the behaviour of the phase centers during 
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a passage above a station as a function of geometric conditions 


- Attitude of the satellite which allows one to relate 
the phase center to the center of mass for the satellite which 
serves as a reference 

- Identical problems for the beacon, which has to be tied 
to a specified material point on the line; the possible change 
of a beacon may not change the tracking network references 

Calibration corrections 


These are Indispensable and fundamental. Thus they have 
to be studied with the same accuracy and have to be available 
for processing: for example. It could happen that a correction 
to be applied will evolve over time. 

Corrections for propagations 

Ionosphere 

The methods of reducing Ionosphere errors are of two 
kinds, If we eliminate the model which appears to still be too 
uncertain : 

- Increase the frequency of operation, because the error 

2 

is proportional to 1/f , 

- operate with two frequencies which have a known ratio /2^ 
so they can be combined to eliminate the Ionosphere effect. 

It Is required that the particular pair selected does not result 
in an amplification of the Instrument errors, if one wishes 
to have a reduction In the Ionosphere error. 

We can assume that it is possible to reduce this error 

_2 

to a level of one centimeter for distance and 10 mm/s for 
speed . 
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Troposphere 


The problem here Is more delicate, and will require a 
detailed In-depth study In order to evaluate the present pos- 
sibilities of correction, because this could have Important 
consequences on the complexity and the cost of the measurement 
apparatus . 

The troposphere error Is Independent of frequency up to 
15 GHz, and Is divided into two components: 

- the dry component which Is modeled relatively well as 
a function of temperature measurements and pressure measure- 
ments at the station (d'Hopfleld model) 

- the humid component which is poorly modeled. For example, 
it would be necessary to measure the water vapor content along 
the wave trajectory, for example, with a radiometer (proposed 

by MacDoran). 

For this troposphere correction, one has to admit that 
one canrot carry out distance and velocity measurements below 
an elevation of about 10 to 15°, or even 20°, If the apparatus 
complexity Is Increased excessively, considering the benefit 
of these low elevation measurements- /25 

2. SELECTION OF THE CALCULATION TECHNIQUE 

Since the errors caused by the Earth potential model are 
assumed to be predominant, we have made an In-depth study of 
these errors. 

EVALUATION OF PERTURBATIONS CAUSED BY THE EARTH'S 

POTENTIAL 


One classical method consists of evaluating the influence 


on osculating elements of terms of the Earth's potential (Kaula 
formalism summarized In the appendix) using a given satellite. 

The existing program has been modified: 

- In order to take Into account the most complete model at 
the present time, GEM lOB (36 degrees, order 36 In the expansion 
of spherical harmonics), 

- In order to calculate the perturbations up to a limit 
of 10 cm. 

It would not serve any purpose to reproduce all these per- 
turbations, but one can summarize the results. 

For a satellite of eccentricity e » 0.001 and 
Inclination of 98^ 5, we obtain - depending on altitude - the 
following number of perturbations: 

at 650 km, 4,100 terms 

at 800 km, 3,100 terms 

atl ,000 km, P,400 terms 

These terms are to be taken into account in the model in 
order to drop down to 10 cm. 

These are results which could be expected, and which 
confirm the problems of a precise orbit calculation. 

Another Important point is the period with which these 
perturbations are produced. 

Approximately, a 5 minute duration corresponds to a 
trajectory length of 2,100 km, and thus one will have sinusoidal 
orbital perturbations with wavelength of 2,000 km, distances 
comparable with that of the oceans. /26 


One finds that the number f'-r which the period is less than 
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15 minutes (6,000 km for the satellite track over the ocean) 
decreases by about 30)5 when the altitude changes from 650 km 
to 1,000 km. 

The greater the altitude, the longer the periods become, 
and In the case of 650 km all of the periods between 5 and 15 
minutes are sampled. 

As an example. Table 1 gives the perturbations In meters 
(for a 17° term and oroer 1) for each of the orbital elements 
with the corresponding period expressed In days or minutes. 

For this term, the summation over p goes from 0 to 17 and the 
summation for q goes from - «» to ♦ « . Only the significant p 
and q values are maintained. Thus one can compare the number 
of terms to be taken Into account for the three different 
altitudes (650 km, 800 and 1,000). 

FUNDAMENTAL CONSEQUENCES 

It seems to us to be almost Impossible at the present to 
have available terrestrial potential models which could translate 
all these perturbations. Thus It Is necessary to reconstruct 
the satellite position In a geometric manner, l.e., to have a 
sufficient variation of measurements In order to localize the 
satellite for a geometric figure. 

Then It is necessary to accumulate these Individual 
locations over time so as to describe the satellite trajectory. 

Using this method, ^>iie Is tempted to Increase substantially 
the number of observation stations of the satellite. There Is 
never the possibility of locating the stations on the oceans. 

For these parts of the trajectories. It Is Impossible to make a 
geometric determination. Thus It is necessary to use satellite 
motion models for describing it. 
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Part of the study thus consisted of evaluating the arc length 
which can be modeled without losing a fixed accuracy. A /28 

certain number of factors play a role, and It Is necessary to 
evaluate them; for example, for networks which observe the 
relative positioning errors, the accuracy, and the measurement 
type, etc. 

GEOMETRIC SOLUTION 

Constraints on simultaneous observation 


- A distance measurement o gives a relation In the form 

p . V (X - ♦ (Y - y)^ ♦ (z - z)^ 

where X, Y, Z are the position coordinates of the satellite 
and X, y, 7 . are the observation station coordinates. 

In order to have the satellite position, three simultaneous 
measurements of distance are necessary for three different 
stations with common visibility of the satellite. 

- In order to obtain the speed of the satellite, we require 
three distinct Doppler effect measurements. In reality, the 
Doppler effect will be translated by a relation of the form 


where and are, respectively, li-e distance at the beginning 
and end of the Integration of the measured Doppler effect. 

Observation relationships 

Overall, for each Instant we will have the following: 

6 observation relationships 

3 distance measurements for the 3 position unknowns 
3 Doppler measurements for the 3 velocity unknowns 
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Obtaining simultaneous tneasurementa /29 

Application of the geometric method requires groups of 
measurements made with 3 beacons or stations, but which correspond 
to a single satellite position. 


In reality, It Is quite probable that the interrogation of 
the beacons will be made In sequence. Thus It Is necessary to 
interpolate among the measurements ci^llected over the entire 
pass above the three beacons . 


Experience has shown that this "smoothing” cannot be made 
with a polynomial with the desired accuracy. It Is preferable 
to select a single motion model for the satellite over the 
duration of the pass. 

The algorithm still has to be defined In order to 
assure that the methodological errors remain smaller than 5 cm, 
and this kind of treatment has to have convergence of the method 
in most cases. /3Q 


INTERPOLATION USING SHORT ARCS 


Principle 

Trajectory observed 
and geometrically 
reconstructed 


Part of trajectory 
connected with 
short arcs 


Trajectory observed 
and geometrically 
reconstructed 
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continent 


We shall assume that for each continent and for the ocean 
coast we have more or less dense observation station networks 
(to be evaluated). These networks will allow the reconstruction 
of the part of the trajectory which Is visible to all 
stations (geometric solution). If we assume that for a given 
arc we can reconstruct che two extremities of the arc under 
consideration, the critical point Is now to evaluate whether 
one can connect them by calculating a trajectory obtained /31 
from satellite motion models, while maintaining the desired 
accuracy. Is It necessary to have an Instrument on an Island? 


The correction between the two extremities will be made 
using the short arc method. 









1. PRINCIPLE 


We shall attempt to determine the perturbations Introduced 
in the orbit when one uses Incorrect dynamic models for /33 
calculating it. This will be done using an ideal measure- 
ment distribution, and then we shall evaluate the Influence of 
a degradation in the orbit coverage by the measurements. 

For this purpose, we shall start by generating a satellite 
orbit as a function of a certain model (often called SIGEDI 
in the illustrations). Prom this orbit, we can simulate 
distance and Doppler measurements for a given collection 
of stations. 

After this, these simulated measurements will be used as 
observations for a trajectory calculation which will be made 
either using a different dynamic model (referenced often by 
the word TRAJEDI in the illustrations), or by introducing 
error models in the station coordinates and/or the measure- 
ments themselves. 

We can divide into two classes the errors which will 
Influence the accuracy with which one can determine the 
satellite position: 

- errors in measurements 

- errors in the force models 

The errors in the measurements Influence the individual 
position of the satellite, whereas the errors in the models 
will introduce errors when one attempts to relate individual 
positions . 

The various steps of the simulation chain and the 
processing chain are described in Tables 3 and ^ respectively. 
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TABLE 3 

SCHEMATIC OP SIMULATION CHAIN 


/ 3 ^ 

Orbit prediction 

- determination of the new simulation period 

- tabulation of the orbit using numerical integration 

creation of a file 

- determination of the collection of visibility periods 
for each station 

- programming of the interrogation Instants 
Generation of meacr^ements , distance and Doppler 

- perfect measurements 

Introduction of noise to the perfect measurements (possibly) 

- real measurements 

Regrouping of the real measurements by station and pass 

Creation of a data bank 

- CHRONOS file (chronological) 

- list of passes 

- data base (measurements) 

Calculation of the station-satellite vector which will be 

used as a reference 


TABLE H 

SCHEMATIC OP DATA PROCESSING CHAIN 


/35 


- Creation of a chronological measurement file and a station 
sub-file 

- Initiation of constants 

- Reading of data 

- Organization of effective parameters (which will be estimated) 

- Calculation of the partial derivative matrix measurements 
with respect to the effective orbital parameters: A 

- Calculation of the contribution of each measurement to the 
second term matrix: B (observed quantity - theoretical quantity). 

- Calculation of the matrix C ■ A ttA (tt * weight) and storing 
of this matrix 

- Inversion of the matrix C and storage of C“^ 

- Possible change of the potential model with respect to the 
simulation 

- Selection of a line of A for each measurement 

- Calculation of the contribution of these measurements to the 

T 

second term matrix D ■ A itb. 

- Solving of the system CX ■ D X ■ C *D, 

- Test of the convergence of the solution in the sense of the 
least squares 

- Adjustment of the effective parameters 

- Calculation of the vector p as a function of the new reconstructed 
orbit 

- Comparison of p and Pq and plot of p - Pq 

2. NUMERICAL VALUES OP THE PARAMETERS /36 

The description of the curves which show the result of 

this study first requires several definitions about the initial 

elements selected. 

For all of the simulations which are used as a reference 
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STATION NETWORK 


/37 


North 

101 

Neskaups tadhur 

latitude T >ngltude 

in degrees 
65,5 - 14 

102 

Jan Mayen 

71.3 

- 8.5 


103 

Scresbysund 

71 

- 22 


104 

Longyearbyen 

78,3 

16 


105 

Hanmer£est 

70,5 

24 


106 

Glasgow 

56 

- 4 


South ’*07 

Ushuaia 

- 55 

. 68 

108 

Rawson 

- 43,5 

- 65 

109 

South Georgia 

- 54,5 

- 37 

no 

Joinville 

- 63,5 

- 55,5 

111 

Montevideo 

- 35 

- 56 

112 

Stanley 

- 57,5 

- 52 


Center 

1 13 

Flores 

- 39,5 

- 

31,3 


1 14 

Canaries 

28 

- 

18 


1 15 

Natal 

- 5,92 

- 

35,16 

Noise 

and 

Bias in the Stations 

and Measurements 



In order to approach the real observation conditions, it 
is possible to introduce random and systematic errors during 
processing. The numerical values used are given below. 


/38 


At 650 and 850 km 
*In the stations: 

. Noise of 20 cm (in this case, we do not desire to Improve 
the positioning of the stations) 
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-5 

. Bias of -1*10 degrees in the stations of the northern 
hemisphere and Natal (in this case, we will attempt to determine 
new coordinates for the group of stations which have a bias) 

At o50 km 


•in the measurements 
. noise of 0.02 m in distance 
. noise of 0.02 mm/s in velocity 

At 850 km 

*In the measurements 

. noise of 0.02 m in distance 

0.02 mm/s in velocity 

. bias of 0.05 m in distance 

-0.1 mm/s in velocity 

*In the stations 

. noise of 0.02m 

-5 

. bias of -1-10 * cos ((^), where 4> is the latitude for the 

southern hemisphere stations and Natal 

General data ,-,r. 


- Simulation step = 60 seconds 

- Interrogation period from one station to another: 
(Interrogation step) = 10 seconds 

- Interrogation period for the same station = 60 seconds 

- Spacing of measurement times: 

- distance date = 4.3 seconds 

- Doppler beginning = 0.4 seconds 

- Doppler end = 4 seconds 

Stated otherwise: 



0.4 4 4,4 8 8,4 12 12,4 
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- Emission frequency (reference) ■ 2 000 000 000 Hz 

- Dlffference frequency ■ 3 750 000 Hz 

- Friction model used is the model derived at 1,000 km, 
1965 with the Barller constants (Jacckia Table). 

Thus the surface to mass ratio is 0,25 . 10*^ o‘/kg. 


Force models 


/^40 


Model GEM lOB 

- Constant used 

Q1 • 0,398600640 D ♦ 15 (m^/s^) 

Earth radius * 0.637814000 D + 07 (m) 

Flattening * 0,298255000 D ♦ 03 

Speed of light = 0,299792500 D ♦ 09 (m/s) 

- Maximum degrees for the terrestrial terms and the 
zonal terms : 36 

Model GEM 8 


- Constants used 

CM - 0,398600800 E 15 
Earth radius = 0,637814500 E + 07 
Flattening = 0.299792500 E + 09 
Speed of light = 0,299792500 E + 09 

- Maximum degrees of the harmonics = 30 

Compromise between GEM lOB and GEM 8 

In order to make the models GEM lOB and 8 compatible, we 
also modified GEM 8 by Introducing the GEM lOB constants given 
above: In the following this Is given as model "GEM 8" with the 
constants of GEM lOB. 


3. NUMERICAL RESULTS 


Description of curves /4l 

The following curves are representative for a selected 
orbit arc; 

R • p * Pq as a function of time 

where p Is he center Earth-satellite vector determined from 

measurements and stations during processing, 

and p Is the vector used as a reference which was calculated 
o 

during simulation, 
p - p Is expressed In centimeters 

The time Is expressed in days, in the Julian calendar. 
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2 . 3 .^ 


SIMULATION OP PERFECT MEASUREMENTS 
TREATMENT WITH THE MODEL GEM lOB 


- 2 . Coordinates of station with noise 
-3 . Coordinates of station with noise 

. Bias In the station coordinates In the South 
and Natal 

. Bias In the station coordinates In the South 
and Natal 


In general, for bias In the coordinates of part of the network 
where the noise over the complex of coordinates has an error 
with a characteristic feature: It Is not the local geometric 
effect which could degrade the Interpretation of altitude 
measurements over a short distance. A noise of 20 centimeters 
Is found almost completely again in P “ Pq (curve 2). 

On the other hand. It seems necessary to have a homogeneous 
precise station network, because a bias over one part of the 
network will lead to a substantial error (curve 3). 

The bias over one part of the network can be reabsorbed, 
which Is indicated and shown in curve 
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/H8 

5.6.7 SIMULATION OF PERFECT MEASUREMENTS 

TREATMENT WITH THE MODEL GEM 10 B BY 
ELIMINATING THE MEASUREMENTS FROM THE TOP 


-5 . Coordinates of stations with noise 

. Bias In the coordinates of the southern stations 
-6 . Bias In the coordinates of the southern station 

-7 . Coordinates of the stations with noise 

. Bias In the coordinates of the southern stations 
(Free coordinates except height) 

If a station Is missing In the center of an arc (Natal In 
this case), we find a slight degradation: 

- comparison between curves 3 and 5» and between curves 
6 and 4 

Nevertheless, It should oe realized that for a single 
pass It Is difficult to want to simultaneously determine the 
satellite orbit and a part of the station coordinates: the 
algorithm to be selected is certainly not the one above, and 
we clearly see the Important role of a preliminary network cal- 
culation. Curve 7 was calculated by not determining the height 
of the stations, and should be compared with curve S where the 
three components were calculated. 
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8-9-10 - SIMULATION OP PERFECT MEASUREMENTS 

- PROCESSING WITH THE MODEL GEM lOB 

- COORDINATES OP STATIONS WITH NOISE 

- BIAS IN THE COORDINATES OF THE SOUTHERN STATIONS 
AND NATAL 

- 9 . MEASUREMENTS OP NATAL ARE ELIMINATED 
-10 . MEASUREMENTS OP NATAL AND CANARY ISLANDS 

ARE ELIMINATED 


The total coverage from the orbit studied Is also an 
Important parameter as can be seen by comparing curves 8, 9 and 
10 . 
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11-12 - SIMULATION OP PERFECT MEASUREMENTS 

- TREATMENT WITH THE MODEL GEM lOB 

- COORDINATES OP STATIONS WITH NOISE 

- BIAS IN THE COORDINATES OP SOUTHERN STATIONS 

- 11 . MEASUREMENTS OP NATAL AND CANARY ISLANDS 

ARE ELIMINATED 

- 12 . MEASUREMENTS OP NATAL AND CANARY ISLANDS 

ARE ELIMINATED AND THE STATION HEIGHT OP 
SOUTHERN STATIONS IS A FIXED PARAMETER 


Another Illustration of the algorithm which can be used. 
Here, we examine the influence of the station altitude during 
processing. 


Distance difference (Earth 
satellite) between 
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- 13 - SIMULATION OP MEASUREMENTS WITH NOISE 
- TREATMENT WITH THE MODEL OEM lOB 


Influence of the measurement accuracy does not produce 
any short period perturbations and has a characteristic 
signature. 
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Distance difference (Earth center/ 
satellite) between 

- orbit tabulated during simulation 

- orbit tabulated in an imperfect case 
i*e« ^ satellite at 6^0 kit 

- *initial parameters of the orbit 

at 7763.37^»305556 

- *Sigedi measurements with noise 

- •Tragedi in the perfect case 




14-15 - SIMULATION OP MEASUREMENT WITH NOISE 

- TREATMENT WITH THE MODEL GEM lOB 

- BIAS IN THE COORDINATE OP THE SOUTHERN STATIONS 
AND NATAL 

- 15 . NOISE IN THE STATIONS COORDINATES 


The mixing of errors in measurements and coordinates 
of stations does not give combinations of the distortions 
which coul d be translated Into short period anomalies. 
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16-17-18 - SIMULATION OP MEASUREMENTS WITH NOISE 

- TREATMENT WITH THE MODEL OEM lOB 

- NOISE IN THE STATION COORDINATES 

- BIAS IN THE SOUTHERN STATION COORDINATES 

-17 . ONLY VELOCITY MEASUREMENTS 
-18 . ONLY DISTANCE MEASUREMENTS 

Here we can see the predominant role of distance measure- 
ments: We find a substantial improvement between curves 17 

(Doppler alone) and 18 (distance alone). One could believe 
that by comparing 16 and 17 that the Improvement due to dis- 
tance is perturbed by the Doppler measurement and, therefore, 
the distance would not be valuable. This Is a concrete Illus- 
tration of the problem of mixing different kinds of measurement 
types where the algorithm to be used must be carefully selected. 
In particular, one has to consider the weighting problems which 
were not into account here. 
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19-20-21-22 - SIMULATION OP MEASUREMENTS WITH NOISE 

- TREATMENT WITH THE MODEL OEM lOB 

- NOISE IN THE STATION COORDINATES 

- BIAS IN THE SOUTHERN STATION COORDINATES 

- 20- . ARC GOES TO THE LEFT; 

VISIBLE ONLY BY 7 STATIONS 

- 21- . ARC GOES TO THE LEFT (7 STATIONS 

visible) 

. ONLY DISTANCE MEASUREMENTS ARE CONSIDERED 

- 22- . ARC GOES TO THE LEFT (7 STATIONS VISIBLE) 

. ONLY VELOCITY MEASUREMENTS ARE 

CONSIDERED 

Influence of the geometry of the path; The path corres- 
ponding to curve 19 is the one shown by a solid line which is 
centered with respect to the stations. The one for curve 20 
is shown by dotted lines on the map. We find a slight degra- 
dation. 

In contast to curves 17 and 18, the predominance of one 
type of measurements with respecc to the other is not clear, 
which Illustrates the geometric factor of these studies very 
well. 
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23-24 - SIMULATION OP MEASUREMENTS WITH NOISE /62 

- NOISE IN THE STATION COORDINATES 

- BIAS IN THE SOUTHERN STATION COORDINATES 

- ARC GOES TO THE RIGHT, VISIBLE ONLY BY 9 
STATIONS 

-23- . TREATMENT WITH THE MODEL GEM lOB 
-24- . TREATMENT WITH THE MODEL GEM 8 


A first estimate of the Influence of the Earth potential 
illustrated by curves 23 and 24. We can Immediately see its 
Important role. The general appearance of the curves up to 
here was regular. The change in the potential model perturbs 
these regular features and shows a kind of undulation around 
the generally regular curves. The evaluation of the very 
short period perturbations (5 to 15 minutes) presages such a 
phenomena which has now been verified. 
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25-26 - SIMULATION OP MEASUREMENTS WITH NOISE 

- NOISE IN THE STATION COORDINATES 

- BIAS IN THE SOUTHERN STATION COORDINATES 

- TREATMENT WITH MODEL OEM 8 

-25- . ARC GOES TO THE RIGHT, VISIBLE ONLY BY 
9 ^STATIONS 

-26- . ARC GOES TO THE LEFT, VISIBLE ONLY BY 
7 STATIONS 


Other illustrations of the influence of potential on 
other arcs. 

The study of the influence of the potential is difficult 
because we have no ac-^ess to reality. The potential models are 
not adequately accurate. 

In order to simulate this Imperfection, we selected the 
change in the model (GEM lOB and GEM 8) which is no less un- 
realistic than to Introduce noise in a given model. All of 
the coefficients are highly correlated in determining the model. 


i 


57 





-27 


- SIMULATION OP PERFECT MEASUREMENTS 

- PROCESSING WITH THE MODEL GEM 8 


/66 


Here, everything Is strictly Identical between the 
simulation and the processing except for a change in the 
potential. Therefore, one can see that it is impossible, 
if one wishes to remain at the 5 cm level, to operate using 
a classical orbiting method. It is therefore necessary to 
at least carry out a partial geometric determination of the 
orbit. 
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Distance difference (Earth center/ | 

satellite between | 

- orbit tabulated during simulation, | 

- orbit tabulated In an Imperfect case, | 

l.e., satellite at 650 km | 

- •initial parameters of the orbit at 7763.37^305556 \ 

- *Tragedl In the perfect case • treatment with GEM 8 i 

I 

'1 

I 
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28-29 - SIMULATION OP PERFECT MEASUREMENTS 

- TREATMENT WITH THE MODEL OEM 8 

-28 , distance measurements only 
-29 . velocity measurements only 


Illustration of the preponderant role of the potential: 
The measurement type Is not significant. 
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30-31 


- SIMULATION OP PERFECT MEASUREMENTS 

- TREATMENT WITH THE MODEL OEM 8 

- ARC COES TO THE RIGHT (9 STATIONS VISIBLE) 

-31- . DISTANCE MEASUREMENTS ONLY 


The Influence of the potential Is not completely 
predictable for two different arcs. 
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32-33 - SIMULATION OP PERFECT MEASUREMENTS 

- TREATMENT WITH THE MODEL GEM 8 

- ARC GOES TO THE LEFT (7 STATIONS VISIBLE) 


-33- . DISTANCE MEASUREMENTS ONLY 
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RESULTS FOR AN ALTITUDE 


Remark: 


OF 850 KM 


Most of the commentary made for the 650 km 
altitude are applicable for 850 km and will 
not be repeated. 
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-3^ - SIMULATION OP MEASUREMENTS WITH NOISE 
- TREATMENT WITH THE MODEL GEM lOB 
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Distance difference (Earth center/satellite) between 

- orbit tabulated during simulation 

- orbit tabulated in an Imperfect case, l.e,, 
satellite at 650 km 

-•initial parameters of the orbit 7760.333333333 
-•Tragedi in the perfect case 





O.31403S 

time in J . 





35-36-37 - SIMULATION OP PERFECT MEASUREMENTS 

- TREATMENT WITH THE MODEL GEM lOB 


-35- . NOISE IN STATION COORDINATES 
-36- . NOISE IN THE STATION COORDINATES 
. BIAS IN THE STATION COORDINATES 
OP THE SOUTH AND NATAL 
. BIAS IN THE STATION COORDINATES 
OP THE SOUTH AND NATAL 
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ematic air in 
of -0.00001 de 
leters of the A 




38-39 


7i4 




- SIMULATION OP PERFECT MEASUREMENTS 

- TREATMENT WITH THE MODEL GEM lOB BY ELIMINATING 
THE NATAL STATION MEASUREMENTS 


-39- . NOISE IN THE STATION COORDINATES 
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40-1*1-^12 - SIMULATION OP PERFECT MEASUREMENTS 

- ELIMINATION OP THE FOLLOW’ STATIONS: 
. NATAL . CANARIES . FLORES 

- TREATMENT WITH THE MODEL GEM lOB 


-40- . NOISE IN THE STATION COORDINATES 
-41- . NOISE IN THE SOU'^H’^D.N STATION COORDINATES 
. BIAS IN THE SOUTHERN STATION COORDINATES 
. BIAS IN THE SOUTHERN STATION COORDINATES 


- 42 - 





4 





-43- . ELIMINATION OP THE NATAL STATION 
MEASUREMENTS 




-*Tragedl in the 
with GEMS 










[ 
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- GTMULATION OF PERFECT MEASUREMENTS 

- ELIMINATION OF STATION MEASUREMENTS: 

NATAL CANARIES FLORES 


- 1 | 5 - 


I 


TREATMENT WITH THE MODEL GEM lOB 
TREATMENT WITH THE MODEL GEM 8 


) 

i' 




*Tragedl in the perfect case 
with GEMS 

*Natal, Canaries and Flores 
eliminated 
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47-48-49 - SIMULATION OP PERFECT MEASUREMENTS 

- TREATMENT WITH THE MODEL GEM 8 

-47- . DISTANCE MEASUREMENTS ONLY 
-48- . VELOCITY MEASUREMENTS ONLY 







00 







<n 








cn 

m 







<n 

o 







tn 








cn 

us 







cn 

4U 








•H 







fO 








CO 



/ 





• 

44 


r 





O 

d 






•» 

sP 

Si 


L 




G 


0 


\ 




G 


44 


\ 




G 

44 

G 


\ 



0 

G 

G 

G 



i 


0 







M 

•H 

44 

44 

44 



\ 

0) 

U 

G 

•H 




\ 

4^ 

to 

G 


4C 




a 


44 

W 





0) 

3 

M 

o 

G 




u 

i 

G • 


G 




*c 

•H 

m 


G 

US 

G 

O 


1 


4J 


•H 

44 





Vi 

to 

O 


44 




(0 

d 

d 

44 

O 




U3 

•H 

G 00 

O 

G 








44 





d 

d 44 

G 

M 




<v 

•o 

•H G 


G 




u 



G 

a 




a 

0 *0 *0 G 

44 





(U 

G G 

G 44 

Si 

G 





G U 

A4 -H 

§ 

X 




0) 

d CO 

G *H 

G 

44 




M-4 

4J iH 

t-4 iH 

U 





M-4 

G 0 

0 G 

G 

d 




•f*4 

JO .0 

UO 44 

a 

•H 




T5 


G G 





mmmimm 


G 4J 

44 G 


•H 




Q> 

AJ 


•0 < 




O 


44 * 

*rl 

G 

2 1 

h- 


c 


•H • 

44 

00 

S 1 

rs 


(0 

iH Jq 

^ G 

•H 

G 





G U 

U • 

d 

U 




0) 

4J O 

O *H 

M 

H 


1 


•H 

G 


<lt 

•K 




Q 

W 1 

1 

1 

1 





\\ 

V 



44 








d 








G 





G G 



B 





G G 



44 





G d 


44 

G 




d 

G G 


G 

G 




o 

44 



U 




•H 

44 G 


44;^. 

44 



0) 

44 

G M 






4J 

G 

G Q 



*K 



d 

iH 

44 






Q> 

3 

V4 


O 

G 



G 

B 

G • 



G 




•H 

G 

tj 


G 

u: 

G 

U 





•H 


44 





00 

O 



44 



to 

d 

d in 


*4-4 

O 



Ed 

•H 

G 00 


o 

G 








44 




0 

d 44 

• 

G 




<V 

*0 

•H G 


U 

G 



O 



iH 

G 

a 



d 

d *0 

•0 <u 

d 

4J 




G 

G G 

G 44 

o 

G cn 

G 



M 

G 4J 

44 «H 


B CO x: 



G 

^ G 

G 

G 

G fO 

44 



44 

G iH 

fH fH 

44 

M CO 




44 

G 0 

0 G 

d 

G fO 

d 00 


•H 

^ 40 

JQ 44 

G 

d-to 

•H 



T3 

^ G 

G G 

e 

CO 




G 44 

44 G 

G 

fH CO 

•H W 


G 

U 


»-i 

G to 

T) CO 


G 

•H 44 

44 

0 

•H 

G 



d 

fH *t4 

•H • 

G 

44 O 

dOUS 


to 

H 40 

Xi G 

G 

•H vP 

G 4 

4 


U 

G }4 

U • 

G 

d r«- 

M -H 

mmm 

CO 

44 O 

O -H 

S *H 

H ; 


C c 

•H 

G 




•K 


9 


G 1 

1 


i 

1 


s 




Q) 

(0 4J 

to *H 
duo 
o o 

1-4 4J f-4 
4J O 0) 
CO 0) > 




50-51 


- SIMULATION OP PERFECT MEASUREMENTS 


-50- . TREATMENT WITH THE MODEL GEM 8 
-51- . TREATMENT WITH THE MODEL GEM 8, BUT THIS MODEL 
TAKES INTO ACCOUNT THE CONSTANTS OP MODEL GEM 
lOB. 


The potential models are determined in a global fashion 
from studying the satellite trajectory, except for GM sometimes 
which comes from planetary probe studies. This Is also true 
for the llllpsold which is used as a reference for the station 
coordinates. One could Imagine that this disparity between GEM 
lOB which was used to simulate the measurements and GEM 8 
with which processing was performed would be the cause of this 
nonregular form (curve 50). Therefore, we gave the same numer- 
ical values as GEM lOB to the fundamental constants of GEM 8 
( curve 51) . 
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52 - 


. SIMULATION OP PERFECT MEASUREMENTS 
. SIMULATION OP PERFECT MEASUREMENTS TAKING 
FRICTION INTO ACCOUNT 

COMPARISON OP TWO RADIUS VECTORS 

-53- . SIMULATION OP PERFECT MEASUREMENTS TAKING 
INTO ACCOUNT 

. TREATMENT WITH THE MODEL GEM lOB, TAKING 
INTO ACCOUNT THE FRICTION MODEL 


In the HASP method, we find that the nongravltational 
accelerations are negligible compared with the effects of the 
Earth's potential. This point was tested with atmospheric 
friction. 

Curve 53 shows the difference between the p calculated 
without atmospheric friction and p calculated with friction. 

The global effect varies between 7 and 9 centimeters over the 
arc length. 

In curve 52, one can see that this effect can be absorbed 
by a coefficient estimated for the past Itself, and, therefore, 
the residual effect is less than 1 centimeter. 
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54-55-56 - SIMULATION OP PERFECT MEASUREMENTS 

- TREATMENT WITH THE MODEL GEM 8 


/ 9 ‘» 


-54- . TABULATION OP THE ORBIT DURING THE FIRST 
TREATMENT ITERATION 

-55- . TABULATION OP THE ORBIT DURING THE FIRST 
TREATMENT ITERATION 

. MODEL GEM 8 TAKES INTO ACCOUNT THE CONSTANTS 
OP MODEL GEM lOB 

-56- . MODEL GEM 8 TAKES INTO ACCOUNT THE CONSTANTS 
OF THE MODEL GEM lOB 


When one tabulates the orbit after the first iteration, 
we can directly see the effect of changing the potential of 
the Earth center-satellite vector (curves 54 or 55). 


These effects cannot be absorbed by an adjustment of the 
Initial orbital conditions: Comparison between curves 55 and 

56. 


88 




57-58-59-60 


-57- . SIMULATION OP PERFECT MEASUREMENTS 
. TREATMENT WITH THE MODEL OEM 8 
-58- . SIMULATION OP PERFECT MEASUREMENTS 
. TREATMENT WITH THE MODEL OEM 8 
. TABULATION OP THE ORBIT AFTER FIRST 
PROCESSING ITERATION 
-59- . SMOOTHING OP CURVE 57 
-60- . RESIDUES BEWEEN CURVES 57 and 59 


Between 57 (after adjustment of tne processing orbit) and 
58 (only effect of changing the Earth potential), we can see tlie 
Inefficiency of the algorithms collected. One could conclude 
that the effects found are from the processing and not from the 
potential . This Is why we smooth curve 57 by curve 59 and the 
differences are plotted In figure 60. This shows the short 
period perturbations due to the Earth potential. 
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61-62-63-64 


-63-62- - SIMULATION OP PERFECT MEASUREMENTS 

- TREATMENT WITH THE MODEL GEM 8, 

MODEL CONSIDERING THE CONSTANTS OP OEM lOB 

-62- - TABULATION OP THE ORBIT AFTER FIRST 

TREATMENT ITERATION 

-63- - SMOOTHING OF CURVE 6l 

- 64- - RESIDUES BETWEEN CURVES 6l and 63 


Results are similar to the preceding case with a simple 
change of the potential. Here, again the short period perturba 
tlons are quite evident In curve 64. 
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65-66-67-68 - SIMULATION OP PERFECT MEASUREMENTS AND 

BIAS IN THESE MEASUREMENTS 

- TREATMENT WITH THE MODEL GEM lOE 

- NOISE IN STATION COORDINATES OP 2 CM 

- BIAS IN THE SOUTHERN STATION COORDINATES OP 
1.10“^ » COS (Oi <l» “ LATITUDE 

-65- . 6 NOR'^HERN STATIONS, 

6 SOUTHERN STATIONS, 

3 CENTER STATIONS 

-66- . i4 NORTHERN STATIONS, 

4 SOUTHERN STATIONS, 

3 CENTER STATIONS 

-67- . 6 NORTHERN STATIONS, 

6 SOUTHERN STATIONS, 

-68- . NORTHERN STATIONS, 

^ SOUTHERN STATIONS 


Here we show the degradation of results when one reduces 
the number of stations (without introducing errors in the Earth 
potential). The suppression of the center otatlons which pro- 
vide complete coverage of the trajectory arc has a substantial 
effect . 


Also note in this case that we introduced a bias in the 
distance measurement and the Doppler measurements. 
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COMPARISON OP THE INFLUENCE OP THE POTENTIAL 
BEIWEN ALTITUDES 650 KM AND 8 50 KM 


We find no substantial difference between the two '•ases, 
except for the general appearance of each curve. 

On the other hand, the type of measurements utilized, 

1 .e . , 

- Doppler and distance 

- distance alone 

- Doppler alone 

does not play any Important role at this level. 




69-70 - SIMULATION OP PERFECT MEASUREMENTS 

- TREATMENT WITH THE MODEL GEM 8 


- 69 - . 650 KM 


Distance difference (Earth center/satelllte) between 

- orbit tabulated durlns slnyilatlon 

- orbit tabulated In an liiq>erfect case, l.e.» satellite at 6S0 km 
-*lnltlal parameters of the orbit 7763.374305556 

■>*Tragedl In the perfect case *treatment with GBM8 


- 69 “ time In & 

Distance difference (Earth center/satellite) between 
“ orbit tabulated during simulation 

“ orbit tabulated in an Imperfect case, i.e., satellite at 850 km 

-♦initial parameters of the orbit 7760.333333333 

“♦Tragedi in the perfect case ^treatment with GEM8 I — iry 


tLue 




71-72 - SIMULATION OP PERFECT MEASUREMENTS 

- TREATMEN. WITH THE MODEL GEM 8 

- ONLY DISTANCE MEASUREMENTS 

-71- . 650 KM 

-72- . 850 KM 





Distance difference (Barth centcr/aatelllte)betveen 

- orbit tabulated during slaulatlon 

- orbit tabulated In an Inperfect case, l.e., satellite at 650 km 
distance measurementa only 

-*lnltlal parameters of the orbit 7763.374305556 i — injn 

~*Tragedl In the perfect case *treatment with GEM 8 | 


- 196 


- 478 



. 759 

“71“ Time In J 

Distance difference (Earth center/satellite) between 
- orbit tabulated during simulation 

• orbit tabulated In an imperfect case, l.e., satellite at 850 km 
distance measurements only 

-♦Initial parameters of the orbit 7760.333333333 

-♦Tragedl In the perfect case *treatment with GEM 8 | 


R 

CN cn. 

- 153 



Time In J 
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73-7M - SIMULATION OP PERFECT MEASUREMENTS 
- TREATMENT WITH THE MODEL GEM 8 


-73- . 650 KM 

-7'»- . 850 KM 
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Distance iilfference (Earth center/aatellita) between 

- orbit tabulated during aimilation 

- orbit tabulated in an iaqi>erfect case, i.e.« satellite at 650 km 
velocity neasurenenta only 

-*initial parameters of the orbit at 7763.374305556 
-’*Tragedi in the perfect case *treatment with GEM 8 f| 
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Distance difference (Earth center/satellite) between 
> orbit tabulated during simulation 

- orbit tabulated in an imperfect case. i.e.. satellite at 850 km 
velocity measurements only 

-^Initial parameters of the orbit at 7760.333333333 I If 

<-*Tragedi in the perfect case ^treatment with GEM 8 | {( 




EM C 

- 156 


- 443 


tlrrie in T*“ 
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/no 

The experience acquired In precise trajectory determination 
allows the formulation of several Ideas for a position determina- 
tion system »t the level of a few centimeters. 

With the reservation of technical feasibility which has not 
been confirmed, but which nevertheless seems possible. It Is 
necessary to define a collection of equipment for achieving 
this object. For example. It is necessary not only to understand 
the Instrument, but also the calculation methods must remain 
compatible with the instrument accuracy. 

The tracking system Includes distance and Doppler measure- 
ments both ways. The measurements are made on board the 
satellite and have only a single time and frequency reference 
for the entire tracking complex. This facilitates the correct- 
ion and transmission of measurements to the computation center. 

This is a concept which allows operations with almost complete 
certainty for any atmospheric conditions. 

Prom the point of view of calculation methods, it is not 
possible to determine an orbit with an accuracy of a few centi- 
meters with a period of several months, using classical techniques. 
This Involves a global reference system for all of the stations 
and an inertial reference system which has to be related to the 
former. This, of course, leads to the development of a method 
based on steadying short arcs (on the order of 1/2 a revolution) 
in which the entire satellite trajectory and the station net- 
work is determined with respect to one another. The station 
network can be defined intrinsically in an Independent manner, 
but will be the only reference with respect to which the orbital 
arcs are calculated. 

/ll] 

Since there is no force model which at the present time 
allows positioning of the satellite with an accuracy of a few 
centimeters and on the other hand, it seems rather unlikely that 
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a global potential model, for example, using ' . cal har- 

monics, could take Into account perturbations or* a few centi- 
meters, we decided to determine the trajectory arc In a geometric 
fashion. This leads to an Increase In the number of observation 
stations: This Is why they are to be designed as simple as 

possible. 

This geometric solution can be simplified because the pre- 
sent analysis shows that it Is not necessary to observe the 
entire satellite trajectory point-by-polnt , but instead, part of 
the trajectory arc can be obsexved either In a partial manner, 
or to Its end. In the latter case, the limit Is on the order of 
2000 to 3000 kilometers. This Is only due tc the short period 
perturbations of the Earth potential because we have shown that 
all of the other error sources (station position, accuracy in 
measurements, nongravltatlonal forces) could not introduce 
short period perturbations. Thus, one can adapt to even the 
ocean geography. 

The concept of the HASP method has been defined and has been 
Justified in broad terms. We believe that its complete analysis 
has not yet been finished. If one really wishes to exploit 
more in depth studies, one will have to determine the arc limits 
to be geometrically observed. Also, algorithms which are best 
suited to this trajectory determination will have to be devel- 
oped , 


If these conditions are satisfied for an ocean or part of 
one, the problems for calculating the complete satellite trajec- 
tory in real time with an accuracy on the order of one or t,<o 
meters will be simply resolved by adding a few stations (on the 
order of 10 at the most) to the network of tracking stations 
covering the region under study. 
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Finally, the simulations made for two satellite altitudes | 

I 650 km and 850 km have not shown any distinct differences, except | 

I that there is a slight degradation for the lowest altitudes when | 

there Is a substantial number of observation stations. We, | 

f therefore, prefer an 850 km orbit but this selection is not | 

I determining for the success of the HASP system. | 
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Development of the perturbation function In osculatlm 


elements . 


We will start with the expression 
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Kaula developed the rollowlng perturbation function 
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Integration of the LaGrange equiitlons 


One Integrates the LaGrange equations by replacing the 
osculating elements s* •* I* ^ In the second terms by 

the average elements a. •. T, !5, a. M where a, e, I are constants 
and Uf w. M are linear functions of time. 


Thus, the total variation of an osculating element is 

equal to the superposition of the variations , to the first 

order produced by each of the terms of perturbation 

function (R - Z . For example, the zonal coefficient C__ 

laicc jO 

gives the following perturbation to the argument w: 


2(l-TVnsHf) '30 




sin 0) 
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In a general manner, the Integration of a term R 

Aa. ■ UK ( ' y . f" * 

n«^''"[a-2p)u «► (L-2?-«s.)M * a(r.-<5)] 


Impq 


1». • U«s 

lagq 
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S, 

where is th;e integral of with respect to its 


argument 



